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SUMMARY 


A theoretical analysis has been conducted to study the 
dynamic behavior and aerodynamic characteristics of a rotor 
in autorotative flight operating in proximity of a reentry 
capsule. A mathematical model representing the rotor entry 
vehicle was established, the necessary equations formulated 
and computer programs were generated for solution of this 
system. 

The equations of motion describing the articulated rotor 
blades retained all nonlinear inertial terms and incorporated 
nonlinear aerodynamics to account for stall, compressibility 
and reversed flow effects. Provisions were also included to 
consider the effects on the rotor of the detached bow shock 
wave generated by the capsule. 

The equations of motion for vehicle dynamic stability 
were written to include coupling effects between capsule body 
motion and rotor forces and moments. Numerical time histories 
of the coupled rotor entry vehicle system are calculated at 
time increments specified by the user. 
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INTRODUCTION 


Background 

Modern aerospace programs have imposed a need for improved 
performance and control capabilities of recovery systems. More 
accurate flight path control during the reentry phase, improved 
maneuverability during descent, lower touchdown velocities, and 
control of the landing mode are sought. 

A concept that has demonstrated capabilities for these 
advanced recovery functions is the unpowered rotary wing de- 
celerator - a multibladed, helicopter-type rotor which provides 
aerodynamic retardation and control in autorotative descent. 
Rotary wing decelerators can be incorporated in a single-stage 
system for low or high speed deployment, retardation and sta- 
bilization, drag force modulation, flight path maneuverability, 
and landing with effectively zero vertical and horizontal speeds. 


State of the Art 

The potential of rotary wing decelerators has been demon- 
strated successfully at subsonic speeds in a USAF-sponsored 
flight test program during which decelerators were dropped at 
speeds up to 550 knots (Reference 1) . A wind tunnel test phase 
of this program also demonstrated that rotors can be deployed 
and operated in axial descent at speeds up to Mach numbers of 
3.0 (Reference 2). Recent wind tunnel investigations by NASA 
(Reference 3) have successfully expanded the operating regimes 
of rotary wing .ecelerators to glide flight at Mach numbers up 
to 3.5. A detailed summary of the state of the art is presented 
in Reference 4. 

The various flight tests and wind tunnel tests mentioned 
above will be enhanced with the complementary use of a rotor 
analysis which will predict blade airloads, blade dynamics, 
rotor performance, and vehicle stability for a rotor entry 
vehicle system throughout its flight spectrum. Available heli- 
copter rotor analyses (e.g.. Reference 5) have been concerned 
with powered rotors in subsonic flight. As such, they are not 
directly usable for studying autorotative rotors in supersonic 
flight; consequently, they must be modified and extended to 
investigate autorot at ing rotors at these extreme operating 
condit ions . 
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Program Goals 


The analysis descrxbed herein was undertaken to px*ovide 
NASA with the analytical tools to evaluate rotor entry vehicle 
systems thereby permitting t-ade-off studies between rotary 
wing decelerators and present day decelerators . Furthermore, 
results of the program will be useful in identifying areas of 
research in rotary wing decelerator technology that require 
further investigation analytically and experimentally. 


Problem Definition 

As mentioned previously, the present state of the art of 
~otor technology must be advanced to examine rotor and vehicle 
dynamic behavior at transonic and supersonic speeds. Theories 
already developed for powered rotors in the subsonic regime 
were used as stepping stones in developing the advanced theories 
and computer programs for predicting the performance, dynamics, 
and stability of autorotating rotors. The key step in the de- 
velopment of the new theories and programs is the formulation of 
a valid blade airloads and rotor performance analysis. Only by 
having reasonable estimates of rotor performance and airloads 
can any confidence be placed in analytical predictions of rotor 
and vehicle system dynamics and stability. 

Consequently, the analyses in the present program are sub- 
divided to evaluate rotor stability, rotor airloads and per- 
formance, and vehicle dynamic stability. Descriptions of these 
analyses and their equations of motion are presented in the 
ensuing sections. 
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ROTOR STABILITY, AIRLOADS, AND PERFORMANCE 


Method of Analysis 

A systematic method of analysis was developed in Reference 5 
for the derivation of the general nonlinear equations of motion 
of a rotor blade with five degrees of blade freedom under both 
steady and maneuvering flight conditions. This method incor- 
porates a position matrix which transforms point mass coordinates 
from a rotating axis system to a fixed axis system by finite 
rotations for each degree of freedom. These position matrices 
are orthogonal. Hinge offsets are obtained by matrix transla- 
t ions . 

In the matrix method of analysis the transformed point mass 
coordinates are differentiated to obtain velocities and accel- 
erations resulting from each degree of freedom. The velocities 
and accelerations are combined to yield the inertial and centrif- 
ugal terms in the blade equations of motion. The blade forces 
and moments associated with the control springs and retention 
springs and damper** are added to these acceleration terms thereby 
resulting in the complete, nonlinear, free equations of motion. 

The forced equations of motion are obtained by applying gener- 
alized forces which incorporate experimental aerodynamic co- 
efficients. 

The aerodynamic characteristics are a function of angle of 
attack, Mach number, and radial station. The unsteady aerodynamic 
forces anu moments are included in the forcing function as an 
approximation. 


Equations of Motion 

Figure 1 shows a deflected, twisted rotor blade assumed *o 
have flapping, feathering, and lagging degrees of freedom. The 
coordinate transformation matrices are obtained by considering 
an axis system initially on the blade feathering axis. The axes 
are then rotated and translated sequentially until they are co- 
incident with the rotor shaft and are translating through space. 
The complete expression for the displacement matrix is obtained 
by combining the coordinate transformations. Velocity and 
acceleration matrices are defined by differentiating the dis- 
placement matrix. The final equations of motion evolving from 
these matrices are derived in detail in Reference 5 and are 
presented below. 
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The inertias and generalized forces in Equations (1) - (3) 
are summarized in Appendix A. For sake of convenience, the 
following symbolism was also used in these equations and will 
be retained throughout the vehicle stability analysis. 

$u=sineu (Du=cosineu 

The feathering input is given by: 


$IN = 0c ^IS 0, c (Cf + 



Numerical Solution 


Numerical methods of solution for the previous differential 
equations of motion are derived in Reference 5. The techniques 
described therein are used to generate response matrices, eval- 
uate transient responses to an initial disturbance, and iterate 
until periodic solutions are acnieved which include nonlinear 
inertial and aerodynamic effects. 

The converged periodic solutions yield azimuthal variations 
of blade kinematics and radial and azimuthal variations of blade 
airloads. These distributions are used ir. the following expres- 
sions to "alcuiate rotor torque, thrust, horizontal and side force. 

Zr R 

Q= ~~ ff H-$*+DG*)rdrd\J/ 

£ - 7r 0 e i 

h 2t r 

T = 2 ^ ff ( LC^+ D$^}(T^drd^ 10) 

o t, 

2 v R 

H / /[H-S^D(D + }$ f -(Ltt*+D$ + )$^0:,]drd^ (6) 


Y = j (-[(Dflq- L$ + ) d, + (L(t,+ D$^) $ s “i,l drdi/z (?) 

°*i 

Because the analysis is written to evaluate an unpowered 
rotor, the equilibrium torque, as defined by Equation (4), must 
be zero. Consequently, at a particular rotor operating con- 
dition, the blade pitch angle must be determined to effect zero 
torque. If the torque does not vanish at a preselected operating 
condition, a torque iteration procedure automatically changes 
blade pitch setting until the torque constraint is satisfied. 


Initially, the method employs a derivative to estimate the 
second trial pitch angle. The value of the derivative is ob- 
tained by using linearized rotor theory (Reference 6) and is 
given by 


dQ 

d e 


7 r 

6 


pQcr\£l R 


■s 


( 8 ) 


Thus, the second trial value of pitch angle is 
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( 9 ) 



Numerical values for the derivative are obtained from the cal- 
culations using the first trial value of blade pitch setting. 


The third trial value of blade pitch angle is estimated by 
using a first-order difference equation_mhich is based on the 


first two calculations. 


Q* - Qo 


l Qz~ Qi \ 

\ 0 1 / 


( 10 ) 


The fourth and subsequent trial values of blade angle are 
determined from the following second order curve. 


Q = a 2 6 Z + a. 6 + a 


(id 


A collocation procedure is used to evaluate the a - coefficients 
in Equation (11). For equilibrium conditions, the torque must 
vanish and new values of Q are obtained from the above quadratic. 


As the torque iteration continues, the three most recent 
values of torque and respective blade angles are retained to 
evaluate the next trial value for blade pitch angle. A toler- 
ance on torque supplied to the computer program as input is 
employed to determine autorotative equilibrium. Retardation 
of the rotor due to bearing friction was omitted from the 
computat ions . 


Supersonic Inflow Through The Rotoi 

A rotary wing that is used as a retardation system on a re- 
entry vehicle will operate through a shock wave generated by the 
capsule during a portion of its trajectory. At transonic and 
supersonic speeds, the rotor blades will be operating in and out 
of the shock wave envelope generated by the capsule. 

For simplicity, the capsule is assumed to be a sphere 
which produces an axially symmetric bow wave as illustrated 
in Figure 2. In order to define which portion of the rotor 
is operating ahead of and which part is operating aft of the 
shock wave, the line of intersection between the detached bow 
shock wave generated by the capsule and the rotor cone must be 
determined . 

The shape of the detached bow wave may be expressed analyt- 
ically using the approximate method of Moeckel as reported in 
Reference 7. Moeckel has indicated that a spherical body 
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generates a shock wave which can be satisfactorily approximated 
by a hyperboloid of revolution with its axis parallel to the 
remote relative wind. The intersection of the hyperboloid and 
the rotor surfaces defines the boundary between disturbed and 
undisturbed flow; the blade flapping response is considered in 
the determination of this boundary. Thus, the line of inter- 
section becomes a function of blade radial station and azimuth 
location. Equations which define the boundary line are presented 
m Appendix B. 

The disturbed flow conditions behind the shock wave are 
evaluated by the oblique shock equations of Reference 8. The 
shock wave angle at rotor azimuth position of 180 degrees is 
used as a control point to determine shock strength. New 
velocities and densities are evaluated behind the shock using 
this control point and the advance ratio and inflow ratio dis- 
tribution are reevaluated. Because of expansion behind the 
shock wave, the air density was also assumed to vary with azimuth 
and rotor angle of attack according to the following relation 

Pf =&,V^(l-.5». R ) + .5[72/>„y4(l + 3(t^ 

-p,V, 2 (l-.5$. n )](t a 2 R (l + C t )}/U 2 (12) 


The effects of flow turning behind the bow wave and the 
effects of mutual shock interference between adjacent blades 
were neglected. 
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VEHICLE DYNAMIC STABILITY 


Equations of Motion 


The six basic equations of motion necessary to describe the 
dynamic behavior of the rotor entry vehicle shown in Figure 3 
are formulated in detail in Reference 9. These basic equations 
are modified to include rotor forces and moments. The modified 
equations are written for a right-hand coordinate system with 
the origin at the vehicle center of gravity. The coordinate 
systems for the vehicle and rotor are shown in Figure 6. All 
rotor hub forces and moments are resolved from the rotor axes to 
the body axes; the rotor axes are defined so that the relative 
free stream velocity lies in the Xr-Zr plane. Furthermore, the 
vehicle was assumed to be symmetrical with respect to the Xr-Zr 
plane thereby eliminating the capsule body cross-products of 
inertia, I xy and I yz . 

The orientation angles between moving earth axes and body 
axes are indicated in Figure 7. The fixed earth axes is also 
portrayed in the figure, and is used as a reference for the 
gravity force. 


A seventh equilibrium equation is required to complete the 
vehicle stability analysis. This latter equation defines the 
rotor angular acceleration in terms of the applied aerodynamic 
torque and the vehicle yaw acceleration. The complete system 
of seven equations is summarized below. 


M B<'Vv B +w Bqe) = F Xf + M B g$e E 4 F Xr 
M glvg-Wgpg+Ugrg) - F, e - + F, R 

M B (wB-u B q B t Vgp B )= F l0 -M B g(I #E C* E +F lR 

PbV'bI.,^ VVb -PB £ lB I « =M ,g +M «R' (B »" B c ) F »r 

-(B M -Bc)F y g 


(13) 

(14) 

(15) 


(16) 



+ A c ( F 2 r + F 2g ) +( Bq) F X g (17) 
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( 18 ) 


^B^-22 * PB^xz + ^yy~ ^xx^ Pb^B + ^xx^b'b ” ^C^ R + ^y^ 

r B -^ = "tT (is) 

^ n 

Methods of Solution 

The preceding differential equations are solved for accel- 
erations in Appendix C. The resulting accelerations are used in 
a self-starting, multistep, predictor-corrector integration pro- 
cedure (Reference 10) to evaluate numerically vehicle velocities 
and displacements relative to the body axes. In turn, these body 
axis kinematic components are integrated to yield translational 
and rotational displacements with respect to a flat earth axis 
system. The equations which describe the earth axis displace- 
ments are developed in References 9 and 11 and are summarized in 
Appendix D. 

In the vehicle stability program, the altitude and free 
stream Mach number are considered constant; thus, the air 
density remains invariant. Furthermore, the entry vehicle is 
assumed axisymmetric and its aerodynamic force and moment char- 
acteristics are assumed specified at a point in the plane of 
symmetry; this point is called the balance center in Figure 3. 
Because the capsule is axisymmetric, body side force, yawing 
moment, and rolling moment about the balance center are neg- 
lected. Consequently, only body lift, drag, and pitching moment 
about the balance center are included in Equations (13) - (18) 
and these aerodynamic characteristics are transferred to the 
body center of gravity. 

In calculating the vehicle stability time histories, the 
rotor airloads Equations (1) - (7) are solved initially until 
equilibrium is achieved in blade airloads, blade motions, and 
rotor torque. These rotor equilibrium hub forces and moments 
are then combined with the vehicle body forces and moments 
according to Equations (13) - (19) and integrated numerically 
at user-specified time increments. After each time increment, 
values of rotor angle of attack, rotor speed, inflow ratio, 
and advance ratio (Q R ,,Q, X , fj. ) are compared to previous 
values. If the increment in these four parameters exceeds 
user-specified tolerances, new rotor airloads are evaluated; 
otherwise, the present rotor airloads axe retained for the 
next time increment. 

Further details of the computer program logic and use are 
presented in the user's manual (Reference 12). 
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CONCLUSIONS 


1. The nonlinear equations of motion are derived which define 
blade transient responses, periodic blade airloads and 
kinematics, and rotor performance of fully articulated 
rotor systems. 

2. Numerical iteration methods are developed to solve these 
equations of motion until equilibrium conditions are 
achieved. 

3. A second numerical iteration method is developed to iterate 
on blade pitch control until the rotor autorotation cri- 
terion is satisfied. 

4. The effect of capsule shock wave interaction with the rotor 
system is estimated by evaluating flow field characteristics 
ahead of and behind the shock wave line of intersection with 
the rotor. 

5. The equations of motion for dynamic stability of the entry 
vehicle are modified to include rotor hub forces and moments 
thereby directly coupling the rotor and vehicle. 

6. All of the preceding analyses are programmed in FORTRAN IV 
language for the Ames Research Center 7040/7094 Direct 
Couple System. 
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APPENDIX A 


ROTOR BLADE INERTIA AND AERODYNAMIC CHARACTERISTICS 


Generalized Inertias 

The generalized inertias used in Equations (1) - (3) are 
defined by the following equations and are grouped into mass, 
first mass moment, and second mass moment integrals. Third 
and fourth order inertial terms are neglected. 


M = f d m 

S, = /(r-ejdm 
S 2 =/£$«* dm 


It =/ (r-e 2 )dm 
I 2 = /(r-e 2 )f S 8x dm 
l3=/(r-e 2 )C(C ex dm 
1 4 =/ £ dm 

^5 = /f 2 ^e x dm 


Generalized Forces 

The generalized forces used in Equations (1) - (3) are 
defined by the following equations. 




where 


<J ^ 

2 \ e o a 
U = 


{[Cl 


+ RU 


dd 


a] (T^+ C 0 $^} dx 


r; U 


an 


(Al) 


(A2) 
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c dC_M ; Co 

RU da C 



dc 

RU da 


a ) (E a + C D S a ]} dx 


(A3) 


Qc 2 _ !^° r 2 c(x-e,) 
I,ft 2 ' 2 c 0 a 




Co (4} dx 


(A4) 


The directions for positive general- 
ized forces are indicated in Figure 4. 

In these equations, it is assumed that the 
lift deficiency function C(k) = 1. The 
local section angle of attack and velocity 
components are shown in Figure £ and are 
given by the following equations. 


a = 0 + + cf> 

U P ~ \(Lp-( x -§2) “jj^ +x{q(Tf-p $*} 

Uj - X2R ^6| + ^(82” ®i ) ( x ~ 62) Cgj ( I * ) ”*■ /f- 

U 2 = U 2 + Up 

* = tan- £ 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
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APPENDIX B 


CAPSULE SHOCK WAVE INTERCEPT WITH ROTOR 


Figure 2 illustrates the rotor entry vehicle at an arbitrary 
angle of attack relative to the free stream velocity. Superposed 
on this vehicle are the approximate sphere and its associated de- 
tached shock wave as described in Reference 7. It is assumed 
that the rotor advance ratio, rotor tip speed, rotor angle of 
attack, and vehicle attitude are known a priori. The free stream 
Mach number can be evaluated from these parameters. 


M 


oo 


/j.£lR 


(Bl) 


If a B = »/2, the free stream Mach number is expressed as 


M 


co 




where C T is the rotor thrust coefficient 


c - I 

" T /t>7TR 2 (nR)2 
and X is the rotor inflow ratio 

\ _ VsinctR-y 

A ~ an 

The hyperbolic shock wave defined in Reference 7 can be 
described in the x 0 , y c, zo wind axis coordinate system as a 
furction of the hyperbola semi-transverse axis (A) and the 
stand-off distance between the sphere and the bow wave { A) • 
Figure 2 shows these dimensions and the appropriate axes; they 
are tabulated below versus free stream Mach number. 


M 

A/Rc 

A/Rc 

1.5 

4.10 

0.540 

2.0 

6.41 

0.340 

2.5 

9.94 

0.225 

3.0 

14.80 

0.180 

3.5 

20.80 

0.160 

4.0 

26.80 

0.150 


At intermediate Mach numbers, A and A can be evaluated by 
linear interpolation. 
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The equation 


in the wind axis system for the hyperboloid is 



(B2) 


In order to describe the line of intersection between this shock 
wave surface and an offset flapping rotor. Equation (B2) must be 
transformed into the rotor axis system. This transformation re- 
quires a standard sequence of translations and rotations as 
described in Reference 5. The resulting expression for the 
hyperboloid in rotor axis coordinates is 


[ X rG«R + (Zr + $aft + ( A + R c + A)] 


[ x RW (z R + B. Q )(tJ+yR 



= 1 


(B3) 


where B e q is the distance between the rotor hub and the center 
of the equivalent capsule sphere. 


The coordinates of a point on the blade are given by 

X R " [ ^2 + ( f ” ^ 2 ) ^ 9 ] 

Yr “ [ e 2 + ( r-e 2^fi] 


z r = ( r ~e 2 )$ fl 


(B4) 


Substitution of coordinates (B4) into Equation (B3) results in 
the following quadratic equation which describes the desired 
line of intersection. 


{cHc!+<gtf t ](M^-l)}(r-e*) 2 

+ 2{C,C 3 -[C 2 C 4 +e 2 (t«$ 2 f ](M^- 1 )} (r-e 2 ) 

+ {Cj-[c5+e|^y(M^,-l )- A 2 } - 0 <B5> 
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where 


C| = + 

Cg = “ So R ^a: f 

Cs= B #q $ aR + ez G„ R (D^+ A + R c + A 

^4 = B 8 q(E aR ~ eg Sa R 


The periodic variation of blade flapping angle with rotor 
azimuth position is known from the numerical solution of 
Equations ( 1 ) - ( 3 ). Thus, the coefficients C^, C2, C3, C4 
and the coefficients of Equation (B 5 ) can be numerically evalu- 
ated at any azimuth position,^/ . The resulting quadratic 
equation at each azimuth can be solved for (r-e2) • 
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APPENDIX C 


VEHICLE STABILITY EQUATIONS 

Equations (13) - (15) can be solved directly for the linear 
accelerations in the body axis system as shown below. 


F x 


U B = 


M e 


B F X r 

— + g $e + Tj + v B r B w bQb 


'e M 


B 


V0= M 


B 


gCe E V M e + W0 P 0 U0rB 


(Cl) 


(C2) 


= M 


Z R 


B 


-g<D #e <t* + j^- + u B q B - v B p B 


8 


(C3) 


Equation (17) can also be solved directly for the body 
pitch acceleration. 


.. M y B + M yR^ B B' B c)F, B + A c (F, R +F r[i ) + (BM-Bc)F, 


< 1 B = 


it 


+ <r|-p|) 


yy 

I„-I, 


« . _ _ ^zi -*-x* 

~ + r B p B ! 


yy 


■yy 


(C4) 


Equation (18) can be written to express yaw acceleration 
in terms of roll acceleration. 


r 8 B - 


Ac( F, R+ F >b ) 


+ <p B -q B r B ) 


xz 


I XX- I 


■zz 


-rz 


xx -Lyy 


•zz 


p B g B 


(C5) 


Substitution of this relation into Equation (16) yields the 
following expression for roll acceleration. 
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I; 


p. 1 6 »- 'v B - - 


-iv>v F .,'} tt: 


Ixx Izz Ixz 
2 


{ 


V j-yy 

iZZ ) izz ixz ^ 

T 

r _ t 2 S 


} 9 B r B + 


Ixz( I XX* I yy Izz ) 


xx i zz ± xz 


Ixx Izz I 


xz 


}p B q 


B 


(C6) 


Finally, the rotor deceleration is simply 

n - r - — 

U ' r B T (C7) 

J-R 

where Q is the torque due to airloads as defined in Equation (4)- 


Because of the axisymmetric characteristics of the capsule 
mentioned m a previous section, the body forces and moments 
given in Equations (Cl) - (C7) are related to the reference 
center lift, drag, and pitching moment by 


Fxg " ( D B a OB Lb 


Fy B = (DeCag-LBSoR) 

F 2B = (D B $ aR+ L B (D BR ) 

M* b = ~ 

My B = My B (Ei, 

Mt B = ^ AcF 7b 

where My B is defined in Equation (C13) . 
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The rotor hub forces and moments in Equations (Cl) - (C7) are 

f„ r = h<c,-y$, 

F, r = Y<t,+H$, 



M >R --(B R -B c ){Y(t, + H$,} + ^ R (t,-M >R $, 

•V (B R -B c )fH(t,-Y$,} + i M (r $, + M, r C, (c 

rv 

where M« and M« are defined by Equations (C16) . 

R y R 

In Equations (C9) and (CIO), the angle tj, between the rel- 
ative wind and the X-axis is defined as follows: 


tj = tan*' 


Vb 

u b 


(Cll) 


The set of equations given by (Cl) - (Cll) can be solved 
numerically after initial conditions are specified and backward 
differences are calculated as described in Reference 10. After 
the pitch and roll rates of the body are calculated from the 
preceding equations, they are transformed into pitch and roll 
rates in the rotor axes according to the following equations. 


q --q B <r,-p B $, 

p = q B VP B c, 

(C12 

The pitch and roll rates from Equation (C12) are then transmitted 
back to the blade equations of motion and used to evaluate new 
equilibrium airloads and blade kinematics. 
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The capsule body lift (Lb) and drag (Dfl) forces are cal- 
culated using the plan area of the capsule body as a reference. 
The body pitching moment (My b ) is also a function of the afore- 
mentioned reference area and capsule body diameter. 


= c l 


2 P 


RC 


Db = 2 P ^ d rc 


My e = 7T/)V« RcC m 


RC 


(C13) 


where the subscript RC refers to the body reference center. 

The rotor hub moments due to the vertical shear forces 
acting at the flapping hinge are calculated in the rotor axis 
system and resolved into the body coordinate system. These 
vertical shears include forces resulting from aerodynamics, 
flapping accelerations and pitch and roll rates. 

The vertical shear acting at the blade flapping hinge at 
any azimuth position can be expressed as: 


dFv 

difr 


+ 2(q$*+ p(t*)] (£0 


(C14) 


The contribution to the moments about the blade flapping 
and feathering axes due to spring and damper effects at any 
azimuth position can be expressed as 






dM, d0 

-K SlN (e-K e ^-K st c-e IN ) 


(C15) 
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The moments due to this force, and these moments resolved into 
the body coordinate system are: 


^!^ =+ 4rr{ (e2 ' e,,$ *' t+ei$ t} + 


d\j/ d \f/ 


dMfi 

d\f/ 


Vf + 


dMe 

d\fr 


(£*- 


+-c 


= ■ -XT’ { (ej ' e| 5 +e ' (t *- t 


d\f/ d yf/ 


, dM* * 

+ ~r,~~ 

df T > (C16) 


» 

Integrating these expressions around the azimuth yields 
the rotor average moments about the body X and Y axes. 


A 

m* r 

A 

M yR 


b ^ 2* dM x 

2 7T J 0 d\jf 

_b r Zw dMy 

2 7T *^o d\j/ 


d\f/ 

d\fr 


(C17) 
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APPENDIX D 


EARTH AXIS DISPLACEMENTS 

The linear and rotational velocities and displacements 
of the rotor entry vehicle in flat earth axis coordinates can 
be calculated from the results of the coupled rotor-body equa- 
tions of motion in Appendix C. The relationships between the 
body velocities and earth axis velocities are given below. 

« 

®E = r B $ » e < D1) 

$£ = P B + (q B V V <te 6 (D2) 

% =( q 8 $ ® E +r B (I ® t J/<I; eE (M) 

Xe = { U B ® BE ' f ( V B^« e ' 1 ' W B ^t^ ^* e 

-{ V 8^ e - W bS* e }S 1 , e (M) 

Ye = { Ub $ep + ( Vb $<*> e + w b (D$ e ) SeJ $* E 

+ {vb^* e " w B $^J (E^ e (d5) 

Ze = ~ 0B$e E +( v b $ 4>e + W gCL^) (L 0e ( D6 ) 


Coordinates in the ea* th axes are calculated by integrating 
Equations (Dl) - (D6! . 
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Nl:r 


ROTOR SHAFT AXIS 



Figure 2. Configuration For Approximating Bow rave Shape 




Figure 3. Vehicle Configuration For Dynamics Study 




Figure 4. 


Generalized Forces Of The Rot*»r System 








Figure 6. Rotor And Vehicle Coordinate Systems 
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